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Abstract 

Uniform estimates for the decay structure of the ro-soliton solution of the 
Korteweg-deVries equation are obtained. The KdV equation, linearized at the 
n-soliton solution is investigated in a class W consisting of sums of travelling 
waves plus an exponentially decaying residual term. An analog of the kernel of 
the time-independent equation is proposed, leading to solvability conditions on the 
inhomogeneous term. Estimates on the inversion of the linearized KdV equation 
at the n-soliton are obtained. 
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1 Introduction 

A little over 100 years ago Korteweg and deVries [O derived their now- 
famous equation 

u t + u xxx - Quu x = 0, 
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in order to explain the observations of solitary waves in barge canals by 
J. S. Russell [[UJ]. Their equation exhibited solitary waves, and so showed 
that such phenomena can plausibly be explained by the Euler equations 
governing the motion of gravity waves in an inviscid nuid.Q This not only 
resolved a long-outstanding controversy, but eventually led to fundamental 
new developments in mathematics. 

In 1967 Gardner, Greene, Kruskal, and Miura showed that the KdV 
equation is in fact a Hamiltonian equation possessing remarkable properties. 
It could be solved by the method of inverse scattering, using the scattering 
Schrodinger equation 

4>xx + k 2 ip - u(x, t)ip = 0. 

As u(x,t) evolves according to the nonlinear KdV equation, the scattering 
data of the associated Schrodinger equation evolves linearly. 

It is still an open problem as to the extent to which solutions of the 
KdV equation approximate solutions of the full Euler equations. There is 
an extensive literature on the validity of the KdV approximation to solitary 
waves of the full Euler equations; cf. for example the works of Amick and 
Kirchgassner [0, Kirchgassner f9j; see also 0, ||. However, in a moving 
reference frame, the solitary wave appears as a time independent solution, 
and is considerably easier to treat mathematically than perturbations of fully 
time dependent solutions of the KdV equation. 

Craig has shown that the KdV equation is a good approximation to 
the full Euler equations over a time period of order e~ 3 , where e 2 = A — 1, A 
being the inverse square of the Froude number (see also Kano and Nishida 
0). It is highly unlikely, however, that the KdV equation is an accurate 
approximation to the Euler equations over an infinite time scale. 

One natural question which one may pose is the following. The KdV 
equation possesses the so-called n— soliton (multi-soliton) solution 



u(x, t) = -2-— r log det ( 5 jk H ) , 9j = ujj(x - atj - 4uA), 

(1.1) 



1 Recently Pego and Weinstein |l4j have discovered that a simple equation equivalent 
to the KdV equation appears in Boussinesq's treatise Q, and, moreover, that Boussinesq 
had noted the existence of solitary waves. 
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where < u% < ■ ■ ■ < u n and acj 6 M, j = 1, . . . , n. 

One may ask, do the multi-soliton solutions of the KdV equation extend 
to full, stable solutions of the Euler equations, say on a semi-infinite time axis 
< t < oo? A positive answer to this question would, as a by product, prove 
the existence of non-trivial time dependent solutions to the Euler equations 
on a semi- infinite time interval < t < oo, something which has not yet 
been done. It would also show that the full Euler equations possess solutions 
which behave like elastically scattering solitary waves. 

The neutral, orbital stability of the multi-soliton solutions of the KdV 
equation has been shown by Maddocks and Sachs [[LJJ based on the observa- 
tion that the n-soliton can be obtained by minimizing the n th conservation 
law of the KdV equation subject to the constraints that the first n — 1 in- 
tegrals of the motion are held fixed. On the other hand, the asymptotic 
stability of the solitary wave has been demonstrated by Pego and Weinstein 
Il4fl , based on a spectral analysis of the linearized Korteweg-deVries equa- 



tion. General arguments, based on the integration of the Korteweg deVries 
equation by the inverse scattering method, imply that the n-soliton solution 
is in some sense asymptotically stable. If so, such a result would be based 
on the analysis of the linearized KdV equation at the n-soliton solution. 

In this paper we analyze the linearized KdV equation at the multi-soliton 
solution. Such an investigation was begun by Sachs |T7[, who constructed a 
representation of the linearized KdV equations using the completeness of the 
squared eigenfunctions of the associated Schrodinger operator. We extend 
his analysis here in the case of the n-soliton solution to obtain estimates in 
norms suitable for studying perturbations of the KdV equation. 

In the case of the solitary wave, one may work in a reference frame moving 
with the wave; the result is that the linearized operator has time independent 
coefficients, and the methods of classical spectral theory of linear operators 
can be applied. In the multi-soliton case this can no longer be carried out, 
and we are forced to consider time dependent operators. Nevertheless, it is 
possible to formulate an analog of the single-soliton follows. 

The n-soliton solution is a 2n parameter family of solutions of the KdV 
equation. Differentiation with respect to those parameters yields a 2n- 
dimensional subspace of solutions of the homogeneous linearized KdV equa- 
tion which decay exponentially in x for fixed time. This subspace plays the 
role of the kernel of the infinitesimal generator in the case of an evolution 
equation with time-independent coefficients, as occurs when the KdV equa- 
tion is linearized about the solitary wave. The presence of this "kernel" leads 
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to 2n solvability conditions. 

We introduce the space W consisting of functions of the form 



n 

ST^ , .2, 



u=2^fj(x-4u;jt) + R(x,t); 
j=i 

the fj are "solitary" -like wave forms which decay exponentially as x — > ±00; 
fj and R(x, t) are analytic in a strip in the complex x plane; and R decays 
exponentially in time as t — > 00, uniformly in x. The class W is closed under 
differentiation and multiplication, an important property when working with 
nonlinear equations. The main purpose of this paper is to invert the KdV 
equation linearized at the n-soliton solution ( |1.1| ) in W. 

In §|2| we prove that the n-soliton solution belongs to W. In fact, for the n- 
solitons, R decays exponentially as \x\ — > 00 as well. This exponential decay 
is not preserved under inversion of the linearized KdV equation. Nevertheless, 
by considering weighted spaces, we can show that the exponential decay as 
x —>■ 00 is preserved. In §[3] we derive the properties of the wave functions 
of the associated Schrodinger operator. In §|4] we give a proof of Sachs' 
completeness theorem, including an extension to a completeness theorem for 
the squared eigenfunctions themselves. 

In §[5] we construct the propagator of the linearized KdV equation. The 
basic estimates on the propagator are obtained in §§|||| In §§ we obtain 
estimates in a Hilbert space of functions analytic in a strip containing the 
real x-axis. In §|8] we obtain estimates in a weighted norm, analogous to 
the estimates in [H| . In §|f] we discuss the inversion of the linearized KdV 
equation in the space W with suitable linear solvability constraints. 

The perturbation scheme of the Euler equations which leads formally 
to the KdV approximation loses derivatives, whereas the inversion of the 
linearized KdV equation is neutral: it is bounded in L 2 , but gains no deriva- 



tives. In a situation such as [jnj the loss of one derivative was compensated 
by the use of weighted norms. Pego and Weinstein used a global existence 
theorem for the generalized KdV equation proved by Kato [Q , together with 
estimates in weighted norms to gain minimal regularity and decay in time. 
The weighted norms exact a toll, however; they do not give rise to Banach 
algebras, and so are difficult to work with in nonlinear problems. 

The loss of derivatives is more severe in the case of the Euler equations, 



and one cannot expect the method in [O] to work. It seems probable that 



some form of hard-implicit function theorem will be needed, such as that 
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described by Moser |T2j or Scheurle JTSJ. A natural space to work in is the 



space of functions analytic in a strip in the complex x plane. 
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2 Asymptotics of the n— soliton solution 

We study in this section the asymptotic properties of the n-soliton solution 
( |1.1| ) of the KdV equation, as t — > oo. The n-soliton solution is a function 
of In parameters Ui, . . . , u n , ati, . . . , a n . Throughout this paper we fix < 
LO\ < • • • < uj n and ax, ... , a n . The speeds of the individual solitary waves 
are 4u \ < • • • < 4u 2 ; the a±, . . . , a n are called the phases. The determinant 
in ( |1.1|) is called the tau function of order n and is denoted by r. 

Theorem 2.1 The n-soliton solution of the KdV equation can be written in 
the form 

n ^2 

u{x, t) = -2j2 ^|sech 2 (^ + 7i ) - 2— log(l + R), (2.1) 



dx 2 



where 



1 1 n ( \ 
In = 77 log(2w n ), jj = - log(2^) + V log ( — — ^ ) , 1 < j < n - 1, 

2 2 k= J+ i 

and 

sup | cosh(ax) R(x,t)\ < Ce~ bt , (2.2) 

x,t>0 

for some a, b > and some positive constant C . 

A similar result is true as t — > — oo, but with different phase shifts 7j. 
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Proof: We will prove that the tau function of order n can be factored 
r(0!, ...,e n ) = 2e-^ + ^ cosh(# n + 7n )r(0 1 + . . . , + /^(l + R n ), 
where r(6 l 1 + /3f , . . . , + f3%_i) * s ^ ne ^ au function of order n — 1, 



0? = log 



! > o, 



and R n satisfies ( |2.2|) . Then by induction, 

n 

r = JJ 2e-^' + ^ ) cosh(% + Ti )(l + #,•), 
i=i 

where each of the Rj satisfies (|2.2j ). The result then follows for the KdV 
solution u upon taking the second logarithmic derivative and letting 



1 + i? = + Rj). 

3=1 



(2.3) 



We begin by writing 

r = det(/ + C n ), C n 



e -{dj+e k ) 

_ 10 j + UJk J i 



i<j,k<n 



As observed in || r can be expanded as a sum of all the principal minors of 
C n . Moreover, each of these principal minors is of the same form. We may 
write 



C n = A n K n A n , A n = diag(e~ ei , . . . , e^"), 

Thus, 
r(0i,...,0 n ) = 



1 + + E e- 2 ^ + ^K^ + --- + e- 2 ^ + - +e ^K^\ (2.4) 



1 



U 3 + U k\ l<j,k<n 



2a;, 

3=1 J l<j<k<n 



where . is the t x t principal minor on the rows ji, ■ ■ ■ ,je- 
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Lemma 2.2 We have 

t(6 u ...,9 n )= T(e u n _O + e-W^T{p x + ft, . . . , # n _! + /£_ x ). 

Proof: In Q2.4| ) the only terms which do not contain the factor e~ 29n are 
the principal minors of the matrix C n _i, so we have the factorization 

n— 1 

r(0 1; . . . , 6 n ) = r(6 u n _ x ) + e" 2 *- (J- + e- 29] K% 

_l_ e -ndj+e k ) K m H \_ e -2(0i+-+e n -i) K (n)\_ 

/ j jkn J 

l<j<k<n-l 



The following formula has been proved in ||, p. 121: 




since 

e -2 7 „ = J_ 

2w„ 
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□ 

From this lemma, we can write the r-function in the form, 9 n = 6 n + j n , 
t(9 u . . . A) =2e^" cosh 6 n ^t{B x , n _i) 

^9 pneti f) 



2 cosh 9 n 

_ r(^ 1 +/?r,...,^i+^_ 1 )) 
2 cosh v„ 7 



e' Un 

+ 



-2e-°« cosh 9 n ( 1+t ^ hgn r(g 1; . . . , 0^) 

+ 1 - t ^ h S (g 1 + ^,...,g B - 1 + /^. 1 )) 
=2e^ coshf n r(0i + ft n , . . . , + + i? n ) 



where 



l-tanh^ l + tanh6> n r(0i, . . . , n _i 

1 + ii„ = 1 



1 + tanh 8„ 
--1 H 



2 r(0! +#,..., n -i + fl£-i) 
r(0i,... A-i) 



To complete the proof of the theorem we have to show that R n satisfies 
Q. Remark first that 



1 < r(0i + /?r, • • • , On-i + PZ-i) < r(9i, 
and that the ratio 

r(0 X) ... A-i) 



r(0i + /3?,... A-i + #U) 

is bounded for all values of 0i, . . . , n _i (it is the ratio of two positive poly- 
nomials of e~ ei , . . . , e -61 ™ -1 having the same degree). 
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Then 



<R ri 



1 + tanh 61 



r^,... A-i) 



r(^ 1 + ^,...,^_ 1 + ^_ 1 ; 



and the required inequality is immediate if x < 0. 
For < x < 2(o£_ 1 + we find; for all a G ^, 

g < <C\e 2Sn = g~ 2tt, nQwg— ax+(2u> n +a)x— Sai^t 

< J^ e -2cu n a ne -ax+[2(wl_ 1 +ul)(2u Jn +a)-8w3]t 

2cu, n 



<_J_ e -2tJ„Q„ e -aa;-[4tj n (^-a;2_ 1 )_2a(£ i ;2+a;2_ :l )]t 



Choose a > 0, a < 2uj n (uj 2 - c^-iV^n + w n-i)> an ^ follows with 

b = 4u n (u 2 - ul_ x ) - 2a{u 2 n + u 2 n _ x ). 

Now, for x > 2{oj 2 1 _ 1 + c<^)t we have 9j = ujj(x — ctj — ^oj 2 t) > —Ujaj, if 
j = 1, . . . , n — 1; and, for < a < 2ui, 

<Rn < r{e x , n _i) - r(0! 0„_i + 



n-l 



: _ e- 2 ^)^- + ^ (1 - e- 2 W+W) e -Wj+^K$ 

3=1 3 l<j<k<n—l 

_| h (1 - e - 2 W + -+l 3 n-l)y-^l + -+9n-l) K (n-l) 



n-l 



n-l 



n-l 



3=1 3=1 



3=1 
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n-1 

aa;-2(w£_ 1 +u£)(2tty-a)t+8w?i 



3=1 
n-1 

<£7 3 e -aa!-[4u'i(^+w^_ 1 -2a;J)-2o(wf +w 2 )]t_ 

i=i 

As before Q follows if < a < 2u;^a; 2 +a;2_ 1 -2u;f)/(u; 2 +u; 2 _ : J. DThis 
theorem shows that the n-soliton solution is asymptotic to a sum of n trav- 
elling solitary waves plus a remainder term that decays exponentially fast to 
zero as t —>■ oo, uniformly in x. 

In the case of 2-solitons, a simple computation shows that 



r = 1 + + + ' - ^ 



2u!i 2uJ2 AUJ1UJ2 \UJ2 + UJi 

We may factor r as 

r = 2e~ (e2+72) cosh(fl 2 + 72 )n 
where 72 = | log(2u; 2 ) and 



-2(01+02) 



n = l + -^e- 2ei (i+tanh(fl 2 + 72 )) + -^ ( U2 Ul ) (l-tanh(0 2 + 72 )). 
As t — > 00, 

(l-tanh(0 2 + 72 )) ->2, (l + tanh(0 2 + 72 )) -f 0; (2.6) 

hence 

n ~ 1 + -L ( ^Z^i) 2 c -a* = 2e~^ coahfa + 7l ), 
2ui \ uj 2 + ooi J 

where 

1 /C(J 2 + 10\ 

71 = - log(2a; 1 j + log 



2 



2 \cg> 2 — cji 

This leads, ultimately, to the factorization 

r = 4e -(0i+0 2 +7i+72) cos h(0 2 + 72 ) cosh(0i + 7i)(l + R), 
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with 



R 



-(1 + tanh(# 2 + 72 ))(1 - tanh(0i + 7l )). 



(2.7) 



Hence, the 2-soliton solution of the KdV equation can be written 

d 2 

u{x, t) = -2uj 2 sech 2 (6 1 + 7l ) - 2cj 2 sech 2 (fl 2 + j 2 ) - 2— log(l + R). 

The diagram below shows the negative of the two-soliton solution ( |1 . 1|) , 
with ui, u 2 = .5, .75, and «i = a 2 = (solid lines). The negative is the lead- 
ing term in the approximation of the free surface in the full Euler equations. 
The first diagram shows the two-soliton solution in the middle of the inter- 
action. The diagram on the right shows the two-soliton solution after the 
interaction. The boxed line shows the superposition of two solitary (sech 2 ) 
waves. These two solitary waves fit the two-soliton solution exactly before 
the interaction. The displacement is the scattering of the solitary waves due 
to the non-linear interaction. 



Theorem 2.3 The n-soliton solution ( |1 . 1|) is analytic in the strip \^sx\ < 
n/2ui n . Moreover, ( |2.1| ) is valid for \^sx\ < n/2uj n and ( |2.2| ) holds for \^sx\ < 
1] , for any < i] < ir/2Q, where Q = max{2 Y^jZi 00 j> ^n}- The constant C 
in Q2.2p depends only on rj . 

Proof: The Gel'fand-Levitan-Marcenko (GLM) equation for the n-soliton 
solutions is (cf. JJ) 





t = 



t = 10 




(2.8) 
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where 



j'=i 

We need to prove its invertibility for complex x. The integral in ( |2.8| ) should 
be understood as a complex integral over r = {z G C : $lz > 9fcc, Qz = Qx}. 
The n-soliton solution is then given by 

u(x, t) = —2—K(x,x,t). 
ax 

Since fl2.8|) is a Fredholm integral equation, its solvability follows from 
uniqueness. We first make a transformation and rewrite the homogeneous 
equation as 

poo 

K x , t (y)+ K Xjt (s)f(s + y + 2x,t)ds = K Xtt (y) = K{x,x + y,t). 
Jo 

(2.9) 

Now it is clear that x appears analytically in the equation, and so the solu- 



tions of (|2.8| ) will be analytic in x wherever we can prove uniqueness. 



To prove uniqueness, multiply ( j2.9| ) by K Xjt (y) and integrate over (0, oo). 
This leads to 

poo poo poo 

0= \K X:t (y)\ 2 dy + / f{s + y + 2x,t)K Xjt {s)K x ~^)dyds 
Jo Jo Jo 

POO n POO 

= \K x j(y)\ 2 dy + y2 e ~ 2UJ]X+8UJ ' t+2UJ3aj \ / ' i2 

Jo j=1 Jo 



e^ s K x As) ds\\ 



The above expression can vanish only when the convex hull of the complex 
numbers e~ 2iu}jV , where x — ^ + irj, contains the negative real axis, hence only 
when 1 77 1 > ir/2uj n . The first part of the theorem is now proved. 

The waveform sech 2 (ujjx) is analytic in the strip \^sx\ < ir/2ujf it has 
poles on the imaginary axis at the points x = i(? + kir)/ujj, fcei Hence 
the left hand side and the sum in the right hand side of ( j2.1|) are analytic in 
the strip < 7i/2u n , so ( |2.1D holds in this strip. 

To obtain the estimate ( |2.2|) remark that the arguments in the proof of 
Theorem ^T] remain valid as long as all the exponential terms e _2 ^ lH 
j = 1, . . . , n— 1 have positive real parts, and l+tanh# n is uniformly bounded. 
These properties hold in any strip \$sx\ < 770, for < 770 < ir/2Q. □ 
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Remark. The argument in the last part of the proof of this theorem im- 
plies also the analyticity of the n-soliton solution in the strip \Qx\ < tc/2Q. 
However, this result is weaker than the one above. We conjecture that fl2.2|) 
holds in fact in any strip \5sx\ < i] Q , for i] Q < Tr/2u n . 

3 Wave functions 

We study in this section some of the properties of the wave functions of the 
Schrodinger operator 

(D 2 + k 2 -u{x,t))ip{x,k) = 0, (3.1) 
when u is an n-soliton solution of the KdV equation 

Ut + u X xx — Quu x = 0. (3-2) 
The Lax pair for this equation is 

L = -D 2 + u, B = -4D 3 + 3(uD + Du). (3.3) 
Substituting u = (p x into the above form of the KdV equation, we obtain 

fx,t + fxxxx — Qfx^Pxx — 0, 

or after integration, the potential KdV equation 

ft + Vxxx - 3ip 2 x = 0. (3.4) 
The linearized KdV equation is 

v t + Vxxx ~ 6{uv) x = 0, (3.5) 
while the linearization of (|3.4| ) is 



ipt + i>xxx - 6uip x = 0. (3.6) 

Thus the linearized potential KdV and the KdV equations are formal adjoints 
of one another. 



Lemma 3.1 Let q be a solution of ( 13.21) , and let fi, f2 be any pair of 
solutions of Lf = k 2 f, f t = Bf, where L and B are given in (|3.3| ). 
Then ip(x,t, k) = (fif 2 )(x,t, k) satisfies the linearized potential KdV equation 



(13.61 ), and ip x satisfies the linearized KdV equation 
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Proof: Straightforward calculation, ^T7\ . □ 

The asymptotic behaviour of the wave functions is determined by the large 
x behaviour of the Lax pair L and B. It is thus given by the simultaneous 
equations 

—D 2 ip = k 2 tp, ip t = -iD 3 ip 

This leads to 

(p + (x, t, k) = m(x, t, k) e - l{kx+ikH \ i{> + (x, t, k) = r{x, t, k)e^ MkH \ 

(3.7) 

The reduced wave functions m(x, t, k), r(x, t, k) are normalized by 

lim m(x,t,k) = l, lim r(x, t, k) = 1. (3-8) 



Lemma 3.2 The reduced wave functions m(x, t, k) andr(x, t, k) are analytic 
in 1 3ft; | < 2tt/Q, and Qk > 0. Moreover, m, r and their derivatives with 
respect to x are uniformly bounded for |3x| < rjo, t G R, Qk > e, for any 
< i]o < 2n/VL and e > 0. 

Proof: The reduced wave functions m(x,t,k), r(x,t,k) are obtained as 
solutions of Volterra integral equations 

/ x I — e 2ik(x-y) 
— u(y, t)m(y, t, k) dy, 
.~ 2ik 



roc j ^—lik(x-y) 

r(x,t,k)=l- I — u{y,t)r(y,t,k)dy. (3.9) 



From Theorem 2.3 it follows that the integrals 



| X _ e 2ik(x-y)\ roc M _ e ~2ik(x-y) I 

\u(y,t)\dy, / — \u{y,t)\dy, 



2\k\ 1 Vtf ' 71 J x 2\k\ 

are uniformly bounded for |3x| < 770, i G M, and SA; > e, for any < 770 < 
27r/f2 and e > 0. Then, the existence, boundedness, and analyticity of m, r 
follow easily from standard arguments for Volterra integral equations. 
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The uniform boundedness of the derivatives of m and r with respect to 
x follows from the uniform boundedness of m and r by using the Volterra 
integral equations (|3.9|) in differentiated form, e.g. 



^-v) u (y,t)m(y,t,k)dy. 



m x (x,t, k) 



□ 



The conjugate wave functions are 



tp~(x, t, k) = (f+(x, t, k), and'0_(a;, t, k) = ip+(x, t, k). 
Since the n-soliton solutions are reflectionless potentials they satisfy 

Lp + = a(k)i/j_, ip+ — a (k)<p_, (3.10) 

where 



a(k) 



W(<p+(x,t,k), j}+(x,t,k)) 
2ik 



For the n-soliton solution of the KdV equation both wave functions 
ip + and ip + are in fact meromorphic functions of k with simple poles at 
—iui, . . . , —iu n . From ( |3.7| ) and ( |3.8| ) we deduce the following principal parts 
expansion of the wave functions: 



ip + (x,t,k) = e-^ kx+4k3t) 



^ + rij(x, t) 
' '•>-■ — ik 



1 U 3 
3=1 



(3.11) 



ij) + (x,t,k) 



i(kx+4k 3 t) 



Pj(x,t) 



3=1 J 



iOn — ik 



(3.12) 



where rij, Pj are the residues of the reduced wave functions at the poles 



k 



-iuj. 



. Then, from ( p.!2j ) and from Lemma |3.2| follows. 



Theorem 3.3 The reduced wave functions m, r and their derivatives with 
respect to x are uniformly bounded in \^Sx\ < i] Q , t G M, $sk > 0, for any 
0<T] < 2ir/tt. 
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We define the eigenfunctions corresponding to the bound states k = iujj 



by 



<fj(x,t) = (p + (x,t,iUj), ipj(x,t) = ip + (x,t,iujj). 



(3.13) 



They can be calculated explicitly by solving the GLM equation, cf. [§]. In 
H the bound state eigenfunctions are normalized differently, namely so that 
ft/)? = 1; this leads to their constraint (2.13). Comparing (2.13) in and 
we see that 



,[GGKM] =d , 



Then, according to (3.3) in ||, the kernel K in the GLM equation is of 
the form 



(3.14) 



K(x,y,t) =-J2e" jaj+ ^%^(x,t>~" jy , 
i=i 

with dj the constants determined above. The substitution of ( |3.14|) into ( |2.8| 
leads to a linear system for tpj\ 



d 3^3 + Yl 



-(6j+9m) 



< x 10 j + U Ti 



-d m ip r , 



(3.15) 



Furthermore, from Theorem 3.4 in || we deduce the following relation be- 
tween the n-soliton solution and the bound state eigenfunctions 



u[x 



t) = -AJ2^e 2 ^^(x,t). 



(3.16) 



We now look for the wave function ip + . For this we determine pj in (|3.12|) 
in terms of ipj. From ( |3.12| ) and fl3.13| ) we find 



ipj (x, t) = e 



—UJjX+iul'H 



Pm(x,t) 



Multiply this equation by dj = e^'" 3 ', and recall that 9j = Uj(x — ctj — 4u?t). 
Comparing the resulting system with ( ft.15 ) we see that 



Pj(x,t) = -e-°idj4>j(x,t) = -e-^ x+4u) P +2u) ^ijj(x,t) 



(3.17) 
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Hence 



n —u)jX+Au)^t+2u)jOij 



cjj — ik 



(3.18) 



We can also obtain an explicit formula for (p + from ( 3.10| ). For this we 
need a(k). First note that 



(p-(x,t,k) =(p+(x,t,k) ~ e l{kx+ik l \ 
Equations ( |3.1Q|) and (|3.18|) then imply 

n 

a{k)= lim e- l(kx+4k3t) ^ + {x,t,k) = 1 - lim V 



as x 



-oo. 



il>j(x,t). (3.19) 



a;,- — zA; 3 
j'=i J 

The limiting values of e~ iipj as x — > — oo can be determined from Q3.15 ). 
Multiply fl3.15| ) by e j and note that e tends to zero exponentially as x — > 
— oo. Moreover, -?/> m is the eigenfunction corresponding to the eigenvalue — 
and is in L 2 (IR). (In fact, it also tends to zero exponentially as x — > ±oo.) 
Hence, these limits satisfy a linear system with constant coefficients, so 



Pn{k) 



10 4 



ik nLi(^i _ ^) 



with P n (k) a polynomial of of degree n. But a(^j) = 0, for j 



,n, 



so 



Pn(£0 =a[J(w j + ife). 

3=1 

The constant a is determined from the asymptotic properties of a(k) for large 
k. From ( |3.19| ) we find 



lim a(k) = 1, 

fc^oo 



so a = (— l) n , and we conclude 



3=1 J 



(3.20) 
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Now, flETUD yields 



3=1 



k — iujj 
k + iujj 



V»i(ar,*) 



(3.21) 



Finally, we compute the coupling coefficients Cj between <pj and ipj, de- 
fined by (pj = Cjipj. From ( |3.18|) a straightforward calculation yields 



n 



CJj — L) m 



2ujj x j- ujj + a; r , 



(3.22) 



Note that, with our choice of normalization of the wave functions, the cou- 
pling coefficients for the bound state wave functions are time independent! 



4 Completeness theorem 



Sachs jjT7|] proved a completeness theorem for the derivatives of the squared 
eigenfunctions of the Schrodinger operator ( |3.1|) , and used this to construct 
the inverse operator for the linearized KdV equation. We describe in this 
section a modification of Sachs' result that gives a direct extension to a 
result for the squared eigenfunctions themselves. For simplicity we assume 
that the potential u is a multi-soliton solution, though the results are still 
valid when u is not a reflectionless potential. 

The squared eigenfunctions satisfy a third order linear homogeneous dif- 
ferential equation, which we write in the form 



Mip = (D 3 - 2(uD + Du) + 4k 2 D)ip(x, k) = 0, 



(4.1) 



where (p is the product of any pair of solutions of ( |3.1[ ). The function u here is 
the n-soliton solution in §||] and §|| (below t will appear only as a parameter, 
and for simplicity we do not write the dependence of the functions on t). 
Using the wave functions (p + (x,k) and ip+(x,k) as a basis of solutions of 
P-lj ) we obtain 



ip 2 + (x, k) , ipl (x, k) , <p>+ (x, k)^) + (x, k) , 
basis of solutions of (Ell). 
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The completeness theorem is based on the following formal calculation of 
the contour integral of the resolvent of the operator M. Let Tr be the semi- 
circle of radius R in the upper half plane traversed from -1 to 1. Formally 

lim — / (D 3 — 2{uD + Du) + 4A; 2 )" 1 8A; dk 

= - lim — / D-\D 2 - 2(u + DuD- 1 ) - A)" 1 d\ (A = -4A; 2 ) 

R-^oo 2-ni J Cr 

=D~\ 

where Cr is the circle of radius 4i? 2 in the A-plane oriented in the counter- 
clockwise direction. 

We construct the inverse of M using as a basis of homogeneous solutions 
the above set of squared eigenfunctions. Since M is skew symmetric, we seek 
a kernel K,(x, y, k) satisfying 

MK(x, y,k) = S(x -y), K(x, y, k) = -K(y, x, k). 

Thus 

[K] = 0, [D X JC] = 0, [D 2 ^] = 1, 

where [/C] denotes the jump of /C across the singularity x = y (from x < y 
to x > y), etc. We also require K,(x,y, k) to be meromorphic in Ssk > and 
bounded as x — > — oo and y — > oo. 
Recall that, for ^sk > 0, 

(p+(x, k) — >0, as x 
if)+(x, k) — >0, as x 
ip + ip + =a(k)i>-i)+ 

We try 

K(x,y,h) = <p%(x, k)Cx(y, k) + (p+(x, k)ip + (x, k)C 2 (y, k), x < y. 

The jump conditions across x = y give three equations for the two unknowns 
C\ and C2, but one finds that the equations are consistent. After some 
computations, we obtain 

K(x, y, k) = ip + (x, k)4>+(y, k) ^S^m ' x<y > 



—00, 
00, 

a(k), as x — > 00. 
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where 

R(x, y, k) = (p + (x, k)i/j + (y, k) - (p + (y, k)ip + (x, k), x <y. 
Based on these considerations, we now prove 
Lemma 4.1 For f G -^i(IR) we have 

8k dk 



(4.3) 



/ f{y)dy= lim / / [K,(x,y,k) - lC(x ,y, 

Jxo R -*°° JTr J-oo 



k)\f{y)dy 



2-ni 



Proof: From the asymptotic properties of the wave functions 
easy to see that 



(4.4) 



it is 



i _ 2ik(y-x) . i 

8kK(x, y, k) = 1 (l + 0(^; 



as \k\ — > oo, 



if x < y. So for X\ < x 2 , 



8k[JC(x 2 , y, k) - K(x u y, k)] ~ < 



( P 2ifc(j/-a; 1 )_ 2ik(y-x 2 ) 

1 xi < x 2 < y, 

e 2ik(x2-y)_ l _ e 2ik(y-x 1 ) _2 



2ik(x2-y)_ p 2ik(x 1 -y) 



xi < y < X2, 
y < Xi < x 2 , 



as k — > oo. For large R, we may replace the integrand in the integral over 
r# by these asymptotic values. Since there is no singularity at k = we may 
deform the contour to an integral over (— R, R) on the real /c-axis. The real 
parts, involving cosines, are odd in k, and hence their contribution vanishes. 
We are therefore left with the identity 

lim - — : / [JC(x2, y, k) — K,(xi, y, k)] 8kdk 

rR sin2k(y — x\) — sm2k(y — x 2 ] 



lim 



R 



2nk 



dk 



[sgn(y - xi) - sgn(y - x 2 )\. 
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The proof of the lemma follows immediately when this result is substituted 
into the right side of fl4~4T ). □ 

The completeness theorem is proved by deforming the integral over Tr in 
(|4.4j ) to the real line. We begin by analyzing the poles of 

8klC = -^ + (x,k)My,k)^^, x<y. (4.5) 
Denote by ' differentiation with respect to k. 

Lemma 4.2 The expression ( |4.5|) has a removable singularity at k = and 
simple poles at k = iuj, j — 1, . . . , n, with residues 

Zujj{a [tujj)) 

for all x, y, where 
Fj(x) = ^ 2 + (x,iUj), 

G j(y) = i( P+(y^ UJ j)-jr( ( P+(y^ k ) - c^+iy.k)) , (4.7) 

UK k=iujj 

with Cj such that (p + (x,iujj) = Cjip + (x,iujj). 

Proof: For k = the wave functions (p + , ip + are real so 

<p-(x, 0) = <p + (x,0), ip-(x,0) = if) + (x,0). 

Then by (|3.10|) we obtain R(x, y, 0) = 0, hence (|4~5| ) is regular at k — 0. 

The poles of 8/c/C coincide with the zeros of a. Hence, by Q3.2C ), it has n 
poles, iujj, j = 1, . . . ,n. From the relation <p+(x, iUj) = Cjip + (x, iuij), where 
Cj is the coupling coefficient defined by fl3.22j ), we see that R(x,y,iujj) = 
and R'(x,y,iujj) ^ 0. Since a 2 (k) has double zeroes at k — iuij, 8k K, has 
simple poles at iuj. A simple calculation shows that the residues there are 
given by 

,^ / x \—MAx,y), x<y , . 

K^x.y) = { y \ yj ' y (4.8) 
\Mj(y,x), y<x 
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where 



Mix v) = V+( X M)1>+(VM) B , X s 



x < y, 



(4.9) 



and 



Bj(x,y) 



dk 



((f+{x, k)tp + (y, k) - ip+(x, k)p+{y, k)) 



x < y. (4.10) 



Using the relationship (f+(x, iuij) = Cjip + (x,iujj), this form of JCj can be 
rewritten in the form ( [4.6| ), as originally obtained by Sachs. □ 

We now deform the contour Tr to the real axis, picking up the residues 
at iujj. By (^4.2| ), the contribution from the real line is 



— / / IC(x,y,k)f(y)dy8kdk 
■<p + (x,k)ip+(x,k) 



ka 2 (k) 



<p + {y,k)i) + {y,k)f(y) dy 



( P+{y,k)ip+(y,k)f(y) dy 



From the equalities ip-(. x >k) — ' l P+( x ,k) = ip + (x,—k) which hold for 
k G E, by Sachs' argument p. 678, this simplifies to 



OO POO 77 

k )^-(y, k) - tf_(x, k)ip 2 + (y, k))f(y) dy ■ ' 

—00 j —00 



Arrik 



(4.11) 



This gives the contribution from the real line; the contribution from the 
residues is straightforward. Hence 



22 



Theorem 4.3 Let u be an n-soliton potential. Then for f G Li(R) ; 

ex n rco 

f(y)dy = Yl / (fcj( x o,y) - fcj(x,y))f(y) dy 

CO j — l J -co 



CO fCO 

l(x,k) - i/j 2 + (x ,k))ilj 2 _(y,k) 

CO J —OO 



dk 

- tfi(x, fc) - ^(x 0l fc))^J(y, fc)]/(y) (4.12) 



where fCj are given in (|4.6|) . 

Theorem 4.4 Lei / and its Fourier transform both be in L\. Then 



co rco 



f(x) = I I [Di> 2 + (x, k)tf.(y, k) - Dtf.(x, ml{y, k)] f(y) dy 



oo J — OO 



Anik 



n rco 

+ J2 C i {DF 3 {x)G J {y)-DG 3 {x)F ] {y))f{y)dy 1 (4.13) 



where Fj, Gj are given in (|4.7| ) and 



Ci = r4 \To - (4.14) 

3 2a; i (a'K0) 2 



Proof: Since /(A:) G Li, / is continuous everywhere, and the derivative 
of the left side of ( |4.12j ) is fix) everywhere. When / and / are in Li, the 
integrals on the right side of ( |4.13| ) converge absolutely. We may divide both 
sides of (|4.12D by (x — xq) and let x — > xq, thereby obtaining (|4.13|) . □ 



5 Propagator of the linearized KdV equation 

Consider the homogeneous initial value problem 

v t + v xxx -6(uv) x = 0, v(x, s) = (f)(x). (5.1) 
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The propagator for this time dependent equation is the operator T t>a such 
that v (t) = T ts cf). By Theorem |44] the solution of the initial value problem 

O is 



v(x, t) 



oo POO 



oo J —oo 



4nik 



[D^l{x,t,k)^ 2 _{y,s,k) 



Di/j 2 _ (x, t, k)ipl(y, s, k)](f)(y) dydk 

n poo 

+ Y< C i [DF j {x,t)G j {y,s)-DG j {x,t)F j {y,s)](j>{y)dy=:T t;s( j>, (5.2) 

-•_ 1 J —oo 



3=1 

where 



Fj{x,t) = ipl(x,t,iuj), 

Gj(y,s) = itp+(y } s,iujj)-^-(ip + (y, s, k) - c^ + (y,s } k)) 

and Cj is given by (|4.14j ). 

A direct computation shows that the functions 

9j(x, t) = ^-(ip + (x, t, k) - Cjip+(x, t, k)) 



satisfy both the Schrodinger equation and the equation ip t = Bip (cf. [[H]])- 
Then by Theorem 3T, Gj satisfies the linearized potential KdV equation 
( |3.6| ) and DGj satisfies the linearized KdV equation ( |3.5| ). The same is true 
for Fj. Hence the last term in ( |5.2|) takes the form of a projection onto the 
kernel of the linearized KdV equation. 

The n-soliton solutions satisfy the KdV equation identically in the In 
parameters «i, . . . ,a n and uj\, . . . , u n . Differentiating the KdV equation with 
respect to each of the a%, . . . , a n and u)\, . . . , u n we obtain In solutions of the 
linearized KdV equation (|3.5|): 



u r 



du 



du 
duij 



j = l,...,n. 



(5.3) 



It is easily seen that the u w . grow linearly with time. 
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We conjecture that the 2n functions DFj and DGj are linear combinations 
of ( |5.3|) . For the DFj, from ( |3.16|) , we have 



Du = -4^2uje 2uJ i a W(rf) = -4^2uje 2 ^ a WFj, 
3=1 3=1 



and by Q 



Hence, 



Tl ^ 

EOU 

J = l J 



y Uj ^DF j = -yp-. 

^ 3 3 A^dai 

3=1 3=1 J 

We conjecture that in fact 

DF = — e~ 2w M — , 
3 Aujj dctj ' 

holds for each j = 1, . . . , n. For n = 1,2 we have confirmed this relationship 
by Maple calculations. 

We expect each of the DGj, j = 1, . . . , n, to be a linear combination of 
all the functions ( |5.3| ). For the one-soliton solution we found, again by Maple 
calculations, 



8^— . flGl = -_ + ^- R j_. (5.4) 

For the two-soliton solution the relationship between DGj and ( |5.3| ) seems 
already to be far more complicated, and we were unable to determine it with 
Maple. 

We see from (|5.4j ) that DG\ grows linearly in t, since the derivative of 
the rz-soliton solution with respect to uj\ has this property. In fact, the same 
is true for all the DGj in the n-soliton be seen by differentiating 

(|3.12| ) with respect to k. This implies that we obtain secular growth terms 
in the sum appearing in ( |5.2| ). In order to eliminate this secular growth in 
the propagator, we must impose the orthogonality conditions 



CO 



G j {y,t)<j)(y)dy= / Fj(y,t)<f>(y) dy = 0, j = l,...,n. (5.5) 



— oo 
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Recall that Fj and Gj are solutions of the the linearized potential KdV 
operator, which is the adjoint of the linearized KdV operator. So, the or- 
thogonality conditions ( |5.5| ) are analogous to solvability conditions in the 
Fredholm alternative for the time dependent operator occurring in (|5.1|) . (cf. 



Lemma \l.2\ below.) 

In the case of a single soliton, one may evaluate the linearized KdV equa- 
tion in a frame moving with the solitary wave. In that case the linearized 
equation is time independent, as are the pair of solutions u ai and u Wl . This 



pair of functions spans the generalized kernel of the linearized operator [TJ 

Lv := v xxx — Aoj\v x — 6(uv) x . 



They are both exponentially decaying as x — > ±00, hence are in L 2 ( 

In the present case, the linearized equation is no longer time independent, 
and we cannot talk about the kernel of a linear time independent operator. 
Nevertheless, the 2n functions ( |5.3|) play a similar role in the analysis. They 
are exponentially decaying in x for fixed time. 



6 Estimates on the propagator 

In this section we estimate the propagator T t a in spaces of functions analytic 
on the strip \$sx\ < 2ir/Q. Given a function <p analytic in this strip, define 
4>a ( x ) — 4>{ x + f° r x ^- Consider the space 

A™ :={0 analytic in \Qx\ < a, ±a G H m (R), 

4> y — * 0, uniformly in \y\ < a, as \x\ — > 00} 

where < a < 27r/fl and m is some positive integer. If a = this space is 
the Sobolev space if m (R). The norm in A™ is defined to be 

\\<f>\\a,m '■= max{|0 a | m , \(j)_ a \ m }, 

where | • \ m denotes the usual norm in H m (M.). 

The following result follows from the properties of the Fourier transform. 

Lemma 6.1 If <p G A™ then the Fourier transform (J) belongs to 

A™ := ($(k) : (1 + \k\) m e a ^(k) G L 2 (R)}. 

Conversely, if <p G A£ then its inverse Fourier transform belongs to A™. 
Moreover, the norms of <ft in A™ and of <p in A™ are equivalent. 
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The norm of G A™ is the L 2 -norm of the function (1 + \k\) m e a ^(f)(k), 
and we shall denote it also by || • || a , m . 
To analyze the propagator we write 

^2 = e a(te+4fc3t)[ 1 + w+ ( a . jfjjfc )] j ^2 = e -x{**+**)[i +w _( x ^ k )]. 

By the results in §[3|, w + is a meromorphic function of fc, with poles in the 
lower half plane at —iuj, j = 1, . . . n, and it decays like 1/k as k tends 
to infinity Moreover, w + is analytic in x in the strip \Qx\ < 2tt/Q] it is 
uniformly bounded in \^sx\ < a, t G R, > — +£, for any < a < 2tt/Q, 
e > 0, and t, fc) = t, k). 

Substituting the expressions for if)± above into the integrand in the first 
term of (15.21), we find 



' [Dij 2 + (x, t, k)^(y, s, k) - Dij 2 _(x, t, k)ij 2 + (y, s, k)} 



Afiik 

= - [cos 2(k(x -y) + 4:k 3 (t - s)) + K^x, y, t, s, k) + K 2 (x,y,t,s, k)] , 

7T 

where 

H^x, y, t, s, k) - H x {y, x, s, t, k) 
Ki{x,y,t,s,k) = , 

H x {x, y, t, s, k) =e 2l{fc( "- y)+4fc3( *- s)) (w + (x, t, k) 

+ w-{y> s, k) + w+(x, t, k)w-(y, s, k)), 

K 2 (x, y, t, s, k) ^[e^-y^^Dw+ix, t, k)(l + w.(y, s, k)) 

- e- 2i ^ x - y)+Ak3{t - s)) Dw_(x, t, £;)(1 + w+(y, s, k))]. 

Since ip+ ( x , t , 0) = ip^(x,t,0), w + (x,t,0) = w_ (x, t,0), and K 2 has a remov- 
able singularity at k = 0. 

The propagator T t)S is thus a sum of 4 terms: 

T = T + T 1 +T 2 + T 3 , 
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where 

poo poo 

^~ J — CO J — CO 



T (t,s)(p=- I I cos 2(k(x-y) + 4k 3 (t- s))<p(y) dy dk, (6.1) 



1 

2^ 



T x (t,s)4>=- / / K x {x,y,t,8,k))<j>{y)dydk, (6.2) 

2 /»oo /* cx: 

T 2 (t,s)4>=- / K 2 {x,y,t,s,k)4>(y)dydk, (6.3) 

^ J —CO J —CO 

n poo 

T 3 (t,s)0=J2 C i {DFjix^Gfaa)- 

j= l J~oo 

DGjfaQFjfasMWdy. (6.4) 
We begin by estimating To. We first note that To can be written 

OO POO 



e 2i(*K*-2/)+« 3 (*-*)) dydk 

— OO J —CO 



1 poo POO 

+ _L / / e-^-v)^*-')) ^y) dydk. 



— CO J —CO 



Changing k — > — k in the second integral, and replacing 2k by k we get 

i /*oo i*co 

T (t,s)<f>=— e Kkx+kHt-s)) e - ik my)dydk. (6.5) 

27T ./ — . / ~v- 



The contribution in the solution of (|7.1|) coming from T is 



/* oo -i poo poo 

Vo (x,t) = - — e ikx e ikt(t-s) e - lky f(y,s)dydkds, 



I ^ J-co 



OO 



or in terms of Fourier transforms 

v (k, t) = - I e ik3 ^ s) f(k, s) ds. (6.6) 
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The action of T (t, s) on is multiplication of the Fourier transform by 
e «fc 3 (t-s)_ This is a unitary operator on A™, so by Lemma |6.1| T is a bounded 
operator on A™, with bound independent of t and s. 

We next show that the operators Tj(t, s), j = 1,2 are smoothing opera- 
tors. The first term in 7\ is in fact a pseudo-differential operator 

T x At, s)<P = — / e 2iikx+4k3{t - s)) w + {x, t, kU(2k) dk. 
2vr J-oo 

From the structure and decay properties in k of w + , Ti i is a bounded map 
from A° a to A l a . In fact, by ( |3.12| ), T^i is a sum of operators such as 



- Pj (x,t) I e^ kx+4k '^-^-^^(2k)dk 

7T k-lUJj 

as well as others which are order k~ 2 . But multiplication of <fi by terms such 
as (k — iuj)' 1 is a smoothing operator; it acts as an integration. Furthermore, 
multiplication by the analytic functions pj(x,t) is also a bounded operation 
on A™. 

We therefore see without difficulty that 

\\Tx,x{t, s)^ <C||0|| a , o , (6.7) 

for some constant C . Thus, by the smoothness of w± in x and y we see that 
Ti t i maps continuously A™ to A™ +1 . The other two terms can be treated in 
exactly the same way. Actually, the third term in T\ is a bounded map from 
A™ to A™ +2 ; though this fact is of no real use. Similarly the operator T 2 is 
a bounded map from A™ to A™ +2 , since its kernel is regular at k = and 
decays like k~ 2 as k — > oo. 

The terms arising from T3 grow linearly in time, as we observed in the 
previous section; and so we have 

Theorem 6.2 Let T ttS be the propagator defined by the initial value problem 
( |5.1| ) ; and let cf) satisfy the orthogonality conditions ( |5.5| ). Then 

I \a,m — C a m j \(f)\ \ a ^m- 
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7 Inversion of the linearized KdV equation 



We now turn to the inhomogeneous equation 

v t + v xxx - G(uv) x = f(x,t), (7.1) 

where u is a multi-soliton solution of the KdV equation. By Duhamel's 
principle, solutions of this inhomogeneous equation are given by 

oo 

v(;t) = ~ J T t>a f(;s) ds. 
t 

The solution of the inhomogeneous equation on the semi infinite interval 
t > is of course not unique. We have chosen this form so that we preserve 
the class of functions decaying exponentially as t — > oo. 
Let 

K™ b = {R(x,t) : supe w ||i2(.,t)lla > m < oo}, 

t>0 

with norm 

||-R|U,b,m = sup e ||-R(-,i)|| a , m , 

where a, m are as in the previous section and b is a positive constant. Denote 
by W the class of all functions f(x,t) of the form 

n 

f(x,t) = Y t f j (x-4uJjt) + R f (x,t), (7.2) 
j'=i 

where Rf e 7Z™i, and fj e ^l™. In fact we assume that fj decays exponen- 
tially to zero as \x\ — * oo, in the strip \$$x\ < a, i.e. fj G A™^, where 

AZ, = {/ : cosh(^)/ G ^}. 

For simplicity we choose the same analyticity domain and the same expo- 
nential decay rate for all fj, though this is not necessary. In this section we 
solve the inhomogeneous linearized KdV equation in the space W. 

The class W corresponds, roughly, to a decomposition of the space of 
functions defined on the half plane t > into exponentially decaying func- 
tions and persisting functions. 
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Recall (Theorem |2.1| ) that the n-soliton solution of the KdV equation is 
of the form 

n 

u(x, t) = Uj(x — 4ujjt) + R u (x, t), (7.3) 
i=i 

where Uj(z) = — 2ujjSech 2 (ujz — ujjOCj + 7j) and R u satisfies Q2.2|). Hence, it 
belongs to W for any < a < n/2fl, < /i < 2u)\, for some b > as in ( |2.2|) , 
and any positive integer m. 

We solve (|7. 1|) in the space W above. For / G W we look for solutions 

i> G W, 

n 

u(;r,i) = ^^(x-4o;|t) + V(x,t), (7.4) 

where V G 7?.™ fe and G ^4™^. Since and Vk decay exponentially in x, 
and are waveforms moving to infinity at different speeds, UjVk G 7?.^ 6 , for 
some b > related to the difference u;| — u;^ of the speeds of and v k] thus 
waveforms moving at different velocities decouple, and their product lies in 
7Z™ b for some b. 

The Vj are determined independently by solving the ordinary differential 
equation 

-4^. + ^ / -6( M ^) , = / r (7.5) 

Then for V we solve 

V t + V xxx -6(uV) x = G(x,t), (7.6) 

where 

n 

G(x, t) = 6 ^(uji;*)* + 6 ^(^)cd + (7.7) 

j^k k=l 

Lemma 7.1 Assume fj G -4.™ M , wzt/i < a < tt/2Q, < fi < 2uj\, and 

/oo 
fj(z)dz = 0. (7.8) 
-oo 
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Then, ( |7.5| ) has a unique solution Vj G A™+ 3 if and only if fj satisfies the 
orthogonality condition 

/•oo 

Uj(z)fj(z)dz = 0. (7.9) 

In particular, if fj(z) = fj(u)jZ — ojjCtj + 7^), and fj is an even function in 
$lz, then ( |7.5| ) has a unique solution Vj(z) = VjiujjZ — UjOij + jj) with Vj an 
even function in ^tz. 



Proof: From (|7.8| ) it follows that there exists gj G A™^ 1 such that fj = 
g'j.f\ Then by integrating ( |7.5|) once, we get 

v" - 4uhj - QujVj = gj. (7.10) 



Recall that Uj(z) = —2u! 2 sech. 2 (ujjZ — ujjdj + 7,). Then, set y 



'jCtj + 7j and (|7.10| ) reads 



UjZ 



] - Avj + 12 sech 2 (y)u,- = h j} (7.11) 



v 

with hj(y) = gj(z)/uj?, hj G A%fy, a' = Uja < vr/2, // = fi/uj < 2. 

In the space above the operator D 2 + 12 sech 2 y — 4 is a compact pertur- 
bation of D 2 — 4, since < \il < 2. It is easily seen that D 2 — 4 possesses 
a bounded inverse from A™, + ^, into A™^~3- Hence the result also applies to 
D 2 — 4 + 12 sech 2 ?/ on the orthogonal complement of its kernel. 

The operator D 2 — 4+12 sech y is well-studied in connection with the 
KdV solitons (cf. fll5p. In L 2 (1R) its discrete spectrum consists of sim- 
ple eigenvalues 5, 0, -3, and its continuous spectrum occupies the interval 
(—00, —4]. Due to translation invariance the derivative of the solitary wave 
is a homogeneous solution of ( |7.11| ), i.e., 



i>o = sech 2 x tanhx, 



is a null function for ( [7.11 ). It spans the kernel of D 2 — 4 + 12 sech 2 ?/ in 
L 2 (R). But 0o belongs to C L 2 (R), so the kernel of D 2 - 4 + 12 sech 2 ?/ 

in this space is also one dimensional. 



2 This is precisely the situation which arises in our perturbation scheme (unpublished 
here) of the Euler equations for surface gravity waves which leads to the KdV approxima- 
tion. 
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From these considerations it follows that ( [7.1 1| ) has a unique solution 
provided hj is orthogonal to 0o- The same is true for the unsealed equation 
( |7.10| ). A simple calculation shows that the orthogonality condition for ( |7.10| ) 
is in fact (|7.9| ) and the first part of the theorem is proved. 

The final part is a consequence of the fact that (TfTTlD can be solved in 
spaces of even functions, and there it always has a unique solution. The 
eigenfunction (j) is odd, so the operator D 2 — 4+12 sech 2 ?/ restricted to the 
subspace of even functions has a trivial kernel. □ 



We substitute now the Vj obtained in this lemma into (|7.7p . Then G G 
1Z™ b1 where b is less than the one given in ( |2.2| ) and 

min{4^(u; 2 +1 - w 2 ), 4//(u;L 1 - a; 2 ), j = 1, . . . , n}. (7.12) 



Lemma 7.2 Assume < a < n /2Vt, and G E TZ™b satisfies the orthogonal- 
ity conditions 

/oo poo 
Gjiy, t)G(y, t) dy = 0, / Fj(y, t)G(y, t)dy = 0, j = 1, . . . , n, 
-oo J — oo 

(7.13) 

where Fj and Gj are the functions in ( |4.7|) . Then ( [7.6| ) has a unique solution 



V(-,t) = -J (T + T 1 +T 2 ) ttS G(-,s)ds, (7.14) 
where Tj, j = 0, 1, 2 are defined in ( |6.1|) -( |OD . 



Proof: Since G satisfies ( |7.13|) , we deduce from the results in §|| that the 
solution V of ( |7.6| ) is given by ( |7. 14| ) . By the estimates in we have 

\\Tj(t, S)G(; s) || a>m < CjWGi; s) || Q , m , j = 0, 1, 2, 

and the constants Cj are independent of t and s. Then 

POO 

\\V\\ a>b ,m =supe 6t / || (T + Ti + T 2 ) t , s G(-,s)\\ ajm ds 

t>0 Jt 



poo 

<Csu P e bt / e~ bs ds\\G\\ aA 
t>o Jt 



C 

in ~^ || G || a,b,my 



33 



hence V £ R™ ft . □ 

The orthogonality conditions (17.13 ) are needed to eliminate terms which 



grow linearly in time. Since Fj and Gj are solutions of the adjoint linearized 
KdV equation, (|7.13j) may be viewed as Fredholm solvability conditions for 



the inversion of the linearized KdV equation on a half space t > 0. 
From these two lemmas we have 

Theorem 7.3 Assume u in ( |1.1|) is an n-soliton solution of the KdV equa- 
tion. Then u is of the form ( [T.3| ), and u £ W ; for any < a < tt/2Q, 
< fi < 2uji, b > as in Q2.2|) , and any positive integer m. 

Assume b is less than the value in ( [7.12[ ) ; and f £ W is of the form ( |7.2| ) 
where each fj satisfies the orthogonality conditions ( [7.81) and ( |7.9|) . T/ien t/ie 
inhomogeneous linearized KdV equation ( |7. 1| ) /ias a solution v £ W ; gwen fry 
(|7.4| ), wrf/i Vj solutions of (\l.5\), provided that 

n 

G(x, t) = 6 ^2(ujV k ) x + 6 + Rf, 

j^k k=l 

satisfies ( |7.13| ). 

In their analysis of the stability of the solitary wave for the generalized 
KdV equation, Pego and Weinstein [jDJ decomposed the equations for the 
perturbation into the domain and range of the linearized operator, and ob- 
tained the orthogonality conditions by imposing modulation equations on 
the phase a and speed u of the solitary wave. Such a perturbation scheme 
is in principle feasible in the case of the multi-solitons; but the fact that the 
perturbation scheme loses derivatives make the problem considerably more 
complicated, as we indicated in the introduction. 



8 Estimates in weighted spaces 

In weighted spaces have been considered to prove the asymptotic stability 



of the solitary wave solutions of the KdV equation. In these spaces the inverse 
of the linearized KdV operator for solitary waves was a smoothing operator. 
In this section we consider similar weighted spaces but for functions which 
are analytic in a strip. We estimate the propagator T and then the solution 
of ( |7.1|) in these spaces. 
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Define the weighted space 

B» v : = {0 : G A™(R)}, (8.1) 
with norm ||0|| a ^ )m := He^^Hc^. Consider also 

B% V fi ■= if(x,t) ■ supe w ||/(-,t)|| a ,„, m < oc}. (8.2) 

t>0 

Denote by || • Ha^m^ the norm in this space. 

From the properties of the Fourier transform one obtains the following 
result. 

Lemma 8.1 If <p G B™ v then its Fourier transform belongs to B™ v := {4>(k) : 

(1 + \k\) m e a W(j)(k + it]) G L 2 (R)}. The spaces B™ v and B™ v are isomorphic 
by the Plancherel theorem. 

The norm in B™ v , being equivalent with the one in B™ v , is also denoted 

by II • II 

To obtain estimates on the propagator T in these spaces we have to 
consider the KdV equation in a moving frame. We shall take the following 
form of the KdV equation 

1 13 
u t - -u xxx + -u x + -uu x = 0. (8.3) 
o 2 2 

This is the actual equation we obtained in the long wave approximation of 
the Euler equations for water waves 

The inhomogeneous linearized equation we solve is 

1 13 

v t - -Vxxx + ~v x + ~(uv) x = f(x, t), (8.4) 

with u an n-soliton solution of (|3.3|) . The results in the previous sections are 
all valid for (gg). 

The solution of (|S.4|) is given by 



0(j 

v(;t) = - J T t)S f(-,s) ds, 
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where T is the propagator of the homogeneous equation. It is defined as for 
(7.1) but with different wave functions ip + and </?+, 

<p+(x,t,k) = e~ i{kx - at) {l + m{x,t,k)), i/)+{x,t,k) = e i(kx - at \l + r{x,t,k)), 
where 

a = -k + -k. 
3 2 

The reduced wave functions m and r are normalized by ( |3.8| ). They have the 
same properties as the reduced wave functions in §|3|. 
The propagator Tq is 



T (t,s) 



1 

7T 
1 

"2^ 



OO POO 



cos 2(k(x — y) — cr(t — s))4>(y) dydk (8.5) 



OO J — OO 



i{ kx+e(s-t)) / e-^^dydk, 



(8.6) 



where a = ^ + |. The other terms in T, Ti, T 2 and T 3 , have the same 
structure as before. 

Choose a such that the decomposition (|2.1| ) and the estimate ( [2.2| ) hold 
in the strip \$sx\ < a. 

Theorem 8.2 Assume < r] < \f?>, b > — 1^(3 — ?? 2 ) ; and m > 1 is an 
integer. Then the propagator Tq satisfies: 



T (t,s)f(-,s) ds 



<a(i7, &)imu 



m— 1,6 



.0- 



(8.7) 



a,r],m,b 



where Ci(rj,b) is a positive constant, and Ci(r),b) — > oo ; as 77 — >• and 
V3. 



Proof: The propagator To from ( |S.5|) is the multiplication of the Fourier 
transform by e la( - s ~^ . Then, to obtain (8.7) we first need bounds for e lCT ( s ~*) 
and fce i<7 ( s ~*), for ^sk = r\. 

Let A; = £ + 277, 0<?7< -\/3, and s > 0. We have 



< e~6' 



and 



36 



where 

' ^ e -^(3-V) s o < s < 4 



/erjs ' jj° 



Hence, we find, for s > t, 

\\T {t,s)f{s)\\ < c q {t],s -t)c l \\f{s)\\ 
for some c% > 0, and by Lemma |8.1| 

\\T (t,s)f(s)\\ < c {v,s -t)c 2 \\f{s)\\ 

a,ri,rn—l- 

Denote 

/oo 
T (t, s)f(s) ds. 

Then 

/oo 
e bt c {r],s - f)||/(s)|| Q)J , jm _ids 

/>oo 

<c 2 / e 6( '- s) c (r7,s -t)ds||/|| aiJ; , m _i i& 



/>oo 

=C 2 / e~ bs C (r], s)ds \\f\\ a , V ,m-l,b = Cl(?7,&)||/||< W n-l,&- 
JO 

It is easily seen that Ci(rj, b) < oo, for < r\ < y/3, and b as in the theorem. 
Furthermore, Ci(f], b) — > oo as rj — > 0, \/3. □ 

The operators Tj, j = 1,2 are again smoothing operators. In fact by 
arguing as in §0 we can prove 



T j (t,s)f(s) ds\\ < Ci(a,77,6)||/|| a , T , )m __ 7 -_i^. (8.8) 
From these results follows. 



Theorem 8.3 Assume rj and b are as in Theorem \8.%\ . If f G B™ b satisfies 
the orthogonality conditions ( |5.5| ) ; then ( |8.4| ) has a solution v G 



2 poo 

W = -E/ T j (t,s)f(s)ds, 

3=0 Jt 
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and 



\ v \\a,T),b,m+l — C||/IU,77,6,-, 



This theorem is the analog of Lemma T71 for the weighted spaces. A 
result similar to the one in Theorem |7.3| holds then for a space W v , defined 



as W but with R f G B£ Vtb , instead of R f G A™ b . 
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